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ABSTRACT 

We present a maximum-likelihood analysis of galaxy-galaxy lensing effects in galaxy 
clusters and in the field. The aim is to determine the accuracy and robustness of con- 
straints that can be obtained on galaxy halo properties in both environments - the 
high density cluster and the low density field. This paper is theoretically motivated, 
therefore, we work exclusively with simulated data (nevertheless defined to match ob- 
servations) to study the accuracy with which input parameters for mass distributions 
for galaxies can be extracted. We model galaxies in the cluster and the field using a 
wide range of mass profiles: the truncated pseudo isothermal elliptical mass distribu- 
tion, the Navarro, Frenk and White profile, and a Power Law model with a core radius. 
We find that independent of the choice of profile the mean mass of galaxies (of the 
order of f 12 M ) can be estimated to within 15% from ground-based data and with 
an error of less than 10% with space observations. Additionally robust constraints 
can be obtained on the mean slope of the mass profile. The two standard parameters 
that characterise galaxy halo models, the central velocity dispersion and the trunca- 
tion radius can also be retrieved reliably from the maximum-likelihood analysis. We 
find that there is an optimal scale i? ma x which marks the boundary between lenses 
that effectively contribute to the measured shear. Lenses beyond -R max in fact dilute 
the shear signal. Furthermore, going beyond the usual formulation, we propose a re- 
parameterisation of the mass models that allows us to put yet stronger constraints on 
the aperture mass of a galaxy halo (with less than 10% error). The gain in signal to 
noise using space observations, expected for instance with the proposed SNAP satellite 
compared to ground based data in terms of accuracy of retrieving input parameters 
is highly significant. 
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1 INTRODUCTION 

Gravitational lensing has now become a popular tool to 
measure the mass distribution of structure in the Universe 
on a range of scales. Recently, there has been consider- 
able progress in mapping the ma ss distribution on rela tively 
large scales using cosmic shear (iRefregier et al.ll2003t) . and 
on cl uster scales combining strong and wea k lensing fea- 
tures (iGavazzi et aljl2003l : iKneib et al.ll2003ft . On the scale 
of individual galaxies as well, there has been much work 
done on modeling and understanding multiple quasar sys- 
tems llFassnacht et al.lll999l : iPhillips et al]l2004l) . In fact, in 
many cases it has become clear that it is almost never a 
unique lens that is responsible for the detected lensing and 
the presence of a nearby galaxy, group or cluster along the 
line of sight plays an important role in inducing the shear 
and amplification llKeeton. Kochanek. fc Seliakll99"^ : lKneibl 



1200(1 iMoller et aljl20o3) . In other words, there are likely no 
clean lines of sight and comprehensive modeling is needed 
to map the lensing configuration accurately. Therefore the 
mass mapping problem is best tackled using an "inverse" 
approach where the adopted method is to model the distri- 
bution of matter around many lines of sight, and optimise 
the mass distribution to match the observations as closely 
as possible. 

Analysing galaxy-galaxy lensing using maximum- 
likelihood lensing techniques is an example of such a method. 
Indeed, the goal of galaxy-galaxy lensing is to obtain con- 
straints on the physical parameters that characterise the 
dark matter halos of galaxies. This is accomplished directly 
using lensing since the deformation in the shapes of back- 
ground galaxies produced by the foreground lenses although 
weak is observationally detected statistically. The difficulty 
is that multiple deflections frequently occur along the line 
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of sight, and therefore, nearby groups or clusters can have 
an important effect yet again on the resultant distorsions. 
This introduces a systematic bias in the mass obtained for 
the deflectors when using simple models. 

Galaxy-galaxy lensing work began with the first de- 
tection of the signal from ground based data (Brainerd, 
Blandford & Smail 1996, hereafter BBS) and later with 
the Hubble Space Telescope (HST) data iGrifntths et alJ 



1996). Maximum- lik elihood techniques have been devel oped 
bvlSchneider fc RbJ (llggjl. iNataraian fc Kneibl l)l99^ and 
iGeieer fe Schneider! il998ti to obtain constraints on galaxy 
halo properties in clusters and in the field. The results 
of these analyses suggest that galaxy halos in clusters are 
significantly less massive but more compact compared to 
galaxy halos around field galaxies of equivalent luminosity 
jNataraianlll998l . lNataraianll2002al) . Besides, in the case of 
galaxy halos in the field no clear edge is detected to the mass 
distribution even on scales of the order of a few hundred kpc 
jMcKav et al.ll200ll;lF7sher et aljEoOpT) . Only two published 
studies to date by Hoekstra (2003); Hoekstra, Yee & Glad- 
ders (2004) have been able to put an upper bound on the 



characteristic extension of a field halo at about 290_ 



kpc and 185l2g h ~ 



h~ 



3^28 kpc, which are only marginally con- 
sistent with each other. Besides these large values do not 
impose a stringent constraint for typical galaxy mass distri- 
butions since at these typical radii the galaxy density is only 
a few times above the mean density of the Universe. 

Galaxy-galaxy lensing studies provide information on 
average properties of the halo population, therefore the re- 
sults depend on the specific parameterised model chosen to 
fit the observational data. From a purely observational point 
of view, the reliability of the galaxy-galaxy lensing signal 
depends on the number density of galaxies whose distorted 
shapes can be reliably measured, as well as any additional 
constraints that can be added to the analysis, for instance, 
redshifts of the lens galaxies, redshifts of the source galax- 
ies, galaxy type, dynamical constraints, and the presence of 
larger scale structure like groups or clusters in the vicinity. 

Other methods to determine the masses of galaxies are 
generally based on the dynamical properties of the lumi- 
nous matter: measurement of the rotation curve or veloc- 
ity dispersion, study of the velocity field of nearby objects 
like planetary nebulae, globular clusters, and satellite galax- 
ies. These dynamical methods are complementary to lensing, 
but often probe much smaller scales. The study of the veloc- 
ity field around galaxies for instance, is generally limited to 
local galaxies, however with the large spectroscopic surveys 
such as 2dF and SPS S, it is now possible to extend such anal- 
yses to larger scales dPrada et al ] |2003l :l Brainerd fc SpeciarJ 
2003). Probing the dynamics of stars in galaxies is limited 
to the inner regions when studying high redshift galaxies. 
Therefore, at the present time there is limited direct overlap 
between lensing and dynamical studies in terms of scales 
probed. Although this situation is likely to change in the 
very near future when large spectros copic surveys of d istant 
galaxies such as the DEEP2 survey foavis et alJl2003l) . the 
VVDS survey ihe Fevre et alJl2003l) or the z-COSMOS sur- 
vey (http://webast.ast.obs-mip.fr/zCosmosl are completed. 

The inner slopes of density profiles provide a strong 
test of structure formation in cold dark matter models and 
lensing provides an unbia sed w ay to estimate the slopes . 
iTreu fc Koopmansl <200-4h and iKoopmans fc Treul <2003h 



have studied the slope of the mass distribution at small radii 
(on scales ranging from a few to about 20 kpc) by combining 
dynamical estimates and strong lensing constraints. They 
find that the mass distribution profile is flatter than the sin- 
gular isothermal sphere (SIS) profile but steeper than the 
NFW profile. Therefore, there is mounting evidence for the 
lack of cores (a constant density region) in galaxies as well 
as in clusters. 

The galaxy-galaxy lensing results from the Sloan Dig- 
ital Sky Survey h ave also provided (iMcKav et alJ 1200 ll : 
ISheldon etai] l2003l interesting constraints on the distribu- 
tion of light and dark matter in galaxies. Mass and light 
trace seem to trace each other reasonably well. The power of 
galaxy-galaxy lensing is that it provides a probe of the grav- 
itational potential of the halos of galaxies out to large radii, 
where no other methods are viable for inter-mediate as well 
as high redshift galaxies independent of the dynamical state 
of the system. A similar approach combining dynamical es- 
timates of the central part of galaxies and galaxy-galaxy 
lensing is planned in the future as part of the GEMS and 
COSMOS projects. 

This paper is organised as follows: in section 2, we de- 
scribe the method adopted to model galaxy lenses, source 
galaxies, the simulations performed to recover the input pa- 
rameters of the lenses and the calculation of the aperture 
mass. In section 3, we present the results for three differ- 
ent classes of lens models considered in this work. In section 
4, we explore the results of re-parameterising the models. 
Whenever necessary our results are scaled to the currently 
preferred flat, low matter density ACDM cosmology with 
Om = 0.3, = 0.7 and a Hubble constant Ho ~ 65 hes 
kms -1 Mpc -1 . In such a cosmology, at z = 0.2, 1" corre- 
sponds to 3.55h^ 5 kpc. 



2 GALAXY-GALAXY LENSING 

We briefly review the basic principles of gravitational lensing 
of distant galaxies before describing the mass distributions 
adopted to model them. This section concludes with the 
presentation of our method to recover the lensing galaxy 
parameters. 



2.1 Lensing equation 

The light rays emitted by a distant galaxy are distorted en- 
route to us by the presence of mass concentrations along the 
line of sight. The distorsion can produce strong effects like 
multiple images or arcs if there is close alignment between 
the distant source and a foreground source, but most of the 
time only a weak distorsion occurs in the galaxy shape. 

Let us consider the multiple lensing equation 
Schneider et alJll992t) . For two lenses A and B, the lens 
equation becomes: 



-Das 
Dos 



Dbs 
Dos 



(1) 



where /3 is the source position, 5a is the deflection due to the 
lens A, <3b is the deflection due to the lens B and Das, D>bs 
and Dos are the angular diameter distances between source 
plane S and lens A, lens B and the observer respectively 
(note that we must have za < zb < Zs). The deflection 
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angle 5x due to the lens X is proportional to the angular 
distance Dox between observer and lens X and to the gra- 
dient of the projected gravitational potential <j> generated by 
the lens X. 

For a given background galaxy (i) and its associated lens 
(j), we can construct the amplification matrix ay, which pro- 
vides the mapping between the source plane and the image 
plane: 



ll 3 
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-7 2 
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(2) 



where k is the convergence, and 71, 72 are the two compo- 
nents of the shear. In the case of multiple deflections (more 
than one lens contributing to the observed distortion), we 
will assume that the total amplification matrix a,: of the 
distant galaxy (i) is equal to the sum of the individual con- 



(3) 



This assumption relies on the fact that we are in the weak 
lensing regime and that the distance between the lenses is 
large compared to the Einstein radius of each individual lens. 
For instance, in the simulations performed in this paper, 
the separation between two lenses is larger than 3" , when a 
typical value for the Einstein radius is about I" . 



2.2 Modeling the mass distribution of galaxies 

Lensing probes the two dimensional projected mass along 
the line of sight, therefore, we deal with the two dimensional 
potential, <f>(R), resulting from the three dimensional density 
distribution p(r) projected onto the lens plane. The related 
projected surface mass density E is then given by: 



47rGE(i?) = V 2 4>(R) 



(4) 



Moreover, we are interested in the two-dimensional pro- 
jected mass inside radius R (the aperture radius R ap er) de- 
fined as follows: 



M aper (R) = 2tt / T,(r)rdr 



(5) 



In this paper, we study three different mass models 
(i) the two component p seudo-isothermal mass distribution 
(PIEMD) iKneibl lT996). which is a more physically moti- 
vated mass profile than the isothermal sphere profile (SIS) 
but sharin g the same profile slope at interm ediate radius; (ii) 
the NFW llNavarro. Frenk fc Whitell997l) profile and (hi) a 
Power Law profile with core radius (PL). These enable us 
to explore a wide range of mass distributions and reveal the 
important parameters that lensing is sensitive to. 



2.2.1 PIEMD profile 

The density distribution for this model is given by: 



p(r) 



Po 



(6) 



(l + r 2 /rI ore )(l + rVr? ut ) 

with the core radius r CO re of the order of 100 pc, and a 
truncation radius r cut . We also introduce a shape parame- 
ter a — r cut /r core- In the centre, p ~ po/(l + r 2 /r 2 ore ) which 



describes a core with central density po. The transition re- 
gion (r core < r < r cu t) is isothermal, with p ~ r~ 2 . In the 



outer parts, the density falls off as p 



, as is usually 



required for models of elliptical galaxies. Fig. illustrates 
this behaviour. These models have been successfully used 
by Natarajan et al. (1998, 2002) to fit observed early-type 
galaxies in cluster lenses. 

Integrating equation |S| we obtain the 2 dimensional 
surface mass density distribution: 



o~ 2 r cu t 



2G(r cut - r core ) y/ri ore + R? y/ r 2 ut + R 2 ' 



(7) 



where gq is the central velocity dispersion for a circular po- 
tential related to po by the following relation: 



1 

^0 /' cut T~ 1 core 

P0 = 2^G ( r 2 r I 

Ail vj / core 1 cut 



(8) 



It can be shown that for a vanishing core radius, the surface 
mass profile obtained above becomes identical to the surface 
mass profile used by BBS for modeling galaxy-galaxy lens- 
ing. The enclosed two dimensional aperture mass interior to 
radius R is: 



M aper (R) 



Tvr cut a 2 . y/r 2 ut + R 2 - Vr'j ore + R 2 s 



(9) 



G f cu t T C ore 

and the total mass of such a model is finite and is given by 



Mtot = 



K<? r cut 



(10) 



G Tcut ~\~ Tcore G 

Fig. shows the behaviour of M aper as a function of 
the aperture radius R ap er and of p(r) for such a profile, 
with (Jo = 220 kms -1 , r core — 0.1", r cut = 30". p(r) is also 
shown, normalised to the critical density of the Universe 
Pcrit, where p crit = 3H 2 /%-kG. 



2.2.2 NFW profile 

The NFW density profile jNavarro. Frenk fc White! Il997ft 
provides the best-fit to the halos that form in N-body sim- 
ulations of collisionless dark matter. In fact the NFW pro- 
file reproduces with good accuracy the radial distribution 
of structures in these simulations over 9 orders of magni- 
tude in mass (from the scale of globular clusters to that of 
massive galaxy clusters). Since it is thought that matter in 
the Universe is dominated by a form of dissipationless cold 
dark matter, this "universal profile" offers an interesting and 
natural way of describing mass concentrations. The density 
distribution of the NFW profile is given by: 



p(r) = 



(r/r s )(l + r/r s 



(11) 



where p s is a characteristic density. It is possible to param- 
eterise this model in terms of M200, which is the mass con- 
tained in a radius r2oo where the criterion p=200p cr it holds, 
and & c the density contrast (or equivalently c = r 2 oo/^"a, the 
concentration parameter). We have the following relations 
between the two parameterisations: 

800 3 „ 200 c 3 .„„. 

p a = d c p c , AI200 = — vrr 2 oopc, o c = — . — r —(12) 

3 3 ln(l + c) - — 



4 Limousin et al. 



The properties of the projected quantities depends on the 
ratio r/r s , so it is useful to introduce the dimensionless ra- 
dial coordinate, x = r/r s = R/r s . Moreover, the velocity 
dispersion <x(r) of this potential, computed with the Jeans 
equation assuming an isotropic velocity distribution, gives 
an unrealistic central velocity dispersion (a(0) =0). In or- 
der to compare the NFW potential with other potentials 
used to model lenses, we define a characteristic velocity a s 
as follows: 



2 ^2 

The surface mass density for the NFW is given by: 



E(x) 
with 

F(x) 



p(r s ,x,z)dz = 2p s r s F{x) 



■\J l — x : 



= argch±) (x<l) 
(x=l) 

= arccos -) (x>l) 



(13) 



(14) 



(15) 



and the two dimensional aperture mass M aDer conta ined 
within the dimensionless radius x is llBartelmanrJll99rj) : 



M apcl (R) = ^iff(x) 



with 



9( x ) 



2G 



In f + ^i=argchi (x<l) 



(x=l) 
(x>l) 




(16) 



(17) 



and the mass M200 can be written as a function of a s , r2oo 
and c: 



M200 = 200tt 



c 2 air 200 



(18) 



S c G 

Fig. shows the behaviour of M ape r as a function of 
the aperture radius R ap er and of p(r), with a s — 225 kms~ 
and r s — 3" . This profile has a concentration parameter c = 
r20o/r a — 12, a typical value for a galaxy, and a projected 
mass inside r2oo: M200 of ~ 3 x 10 12 Mq . 

2.2.3 Power Law profile with a core 

Another simple model to describe the mass distribution of 
a galaxy is a power-law model (PL) with a core. In a CDM 
dominated hierarchical structure formation scenario, mass 
profiles are expected to be independent of the mass scale, 
therefore a power-law profile is of interest. The PL mass dis- 
tribution has three parameters: a core radius r core of the or- 
der a kpc for an average galaxy, a central velocity dispersion, 
00 measured in (kms -1 ) and an exponent (a) which defines 
the gradient of the mass distribution. The three-dimensional 
density profile is (Kneib 1993): 



P(r) = po 



(! + (7-37) 2 ) 2+c 



(19) 



Introducing x — — - — , the density profile falls off as p ~ 

x~ 2 ^ 1+a \ Note that the case a=0 corresponds to an isother- 
mal sphere with a core radius, and a >0 defines density pro- 
files steeper than an isothermal sphere with a core radius. 
The relation between po and ao is given by: 



al 9(1 -2a) 

po ~ V 2 d^r — 

' core v 1 

The surface mass density is: 



1 + ^x 2 
(l + x 2 ) 3 / 2 +° 



(20) 



(21) 



and the two dimensional aperture mass contained within the 
dimensionless radius x is: 



Maper (R) 

with 

Il+a 



3(1 + 2a)alr core x 2 + x A 



G 



(! + 2-2)3/2+, 



--Ti+o 



da; 



(22) 



(23) 



(l + x 2 y+ a 

Fig- shows the behaviour of p(r) and M aper {r), for 
such a profile, with a = 400 kms~ and a = 0.4 (in this 
case, 7i.4=1.06). 

2.3 Comparing these profiles 

Beginning with the Boltzmann equation to describe the be- 
haviour of the cold dark matter collisionless particles that 
constitute a galaxy halo, the Jeans Equation for a spherical 
potential and for an isotropic velocity distribution (a = a r ) 
is given by: 

1 d{a 2 p) _ _d£ 
p dr dr 



(24) 



where $ is the three dimensional potential. Considering the 
limit at +00, wherein a(r) — > and p(r) — > to perform 
the integration, we have: 



2, s 1 r 00 1 > ^( r ') , > 



(25) 



Solving this equation, does not lead generally to a simple 
analytical expression for the velocity dispersion. 

Hence, one usually uses the rotation velocity defined as: 

GM aper {R) 



R 



(26) 



Fig-Qshows the behaviour of V ro t{R) for the three dif- 
ferent mass profiles studied here. 

It can be shown that for any spherically symmetric pro- 
file, the mass inside the Einstein radius Re is proportional 
to R 2 E , so that profiles for which Re is constant can be com- 
pared. It is easy to show {e.g. Kneib 1993) that: 



M aper (R E 



■ 7rE cr it Re 



(27) 



Thus we adjust the parameters of the different mass pro- 
files in order to have the same Einstein radius, and therefore 
the same mass within the Einstein radius. The results are 
illustrated in Fig. [5] this plot allows us to rescale the veloc- 
ity dispersions derived for each profile. We can immediately 
see from this plot that for r > 4" , the shears computed from 
the three profiles are comparable. Note that in the case of 
the PL profile there is a strong dependence of the velocity 
dispersion on the value of the exponent a. For example, the 
PL profile with ao = 660 kms -1 and an exponent of 0.3 in- 
duces the same shear as one with ao = 400 kms -1 and an 
exponent of 0.4. Therefore, there is degeneracy between the 
value of ao and a for the PL profile. In order to illustrate 
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Figure 1. The density profile p(r), the rotation velocity V ro t(R) 
and the aperture mass M a per(R) for the 3 different mass profiles 
studied in this work. For each of these model profiles the relevant 
parameter choices are shown in the figure. 
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Figure 2. The reduced shears (j-z— ) for the 3 different models 
for which the Einstein radius Re = 1" (upper panel), and dif- 
ferences between them expressed as a percentage where we have 
used the PIEMD as the reference profile (lower panel). 



the behaviour of the PL profile, we include this latter profile 
in Fig. Hand |U 



3 SIMULATING GALAXY-GALAXY LENSING 

3.1 Measurement of background galaxy shapes 

We study in detail the lensing effects in two observational 
scenarios: (i) a ground based survey using a wide field cam- 
era and (ii) data from space-based observations. Ground 
based data is characterised by the following image quality: 
seeing of about 0.6-0.8 arcsec, and a galaxy number density 
of 20-40 galaxies per square arcminute, of which only 50- 
70% can generally be used to measure their shapes reliably. 
The above estimates were obtained from two hours of obser- 
vation in the R band with the CFH12k camera with a field 
of view of 44x28 square arcminutes. 

Space observations have a significantly better image 
quality, with a PSF of about 0.1 to 0.15 arcsec, and a 
galaxy number density of about 40-100 galax ies per square 
arcmi nute (SNAP mission sensitivities - see iRhodes et alJ 
( 2003)), of which about 80%. can be used in a weak lensing 
study since their shapes can be measured to the requisite de- 
gree of precision. Half an hour of observation in the R band 
of A2218 with the HST have provided the above estimates. 

The measurement of shapes of lensed background galax- 
ies is made using the second moment of the intensity of their 
light distribution. The quantity that is extracted for each 
galaxy is its complex ellipticity, e defined as e = e\ + ie-z. 
The magnitude of the ellipticity is e = yef + ef) the 



mean ellipticity m 



£ 4 < 



the dispersion of the elliptic- 
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ity is equal to the square root of the variance, defined by 
var = ' N , where N is the number of objects. 

3.2 Scaling the mass distribution from the light 
distribution 

The foreground lenses are described by a mass profile with 
known input parameters. The parameters used to describe 
the different lenses are scaled as a function of luminosity. 
The scaling relation for ero assumes that mass traces light, 
and its origin resides in the Tully-Fisher or Faber- Jackson 
relations. The scaling relation for the radial parameter as- 
sumes that the mass-to- light ratio is constant for all galaxies. 
Note that they are other possible scaling relations, and that 
in principle we can test them with lensing. 

3.2.1 PIEMD profile 
We have for this profile: 

ff o = <xS(-^)* & r cut = r* ut (-^)$ (28) 

The parameter r core is kept fixed at 0.1" - a fairly typical 
value for a galaxy. 

From equation (13), we can scale the total mass with 
the luminosity as: 

2 *2 * t 

. , 7rcr r et ,t na r cut L , 3 / 4 

Mtot - q - q (-j— ) (Jy) 

3.2.2 NFW profile 

Similar to the PIEMD profile, we have: 

°o=Va(jj;) i & r s =rt(^ (30) 

3.2.3 PL profile 
For the PL profile: 

°"o = o{j^) T & r core = r* ore (-^)' (31) 

In order to illustrate the coherence of these scaling laws, 
we show the shear profiles obtained for a typical faint (L = 
L * /10) and bright (L — 3L*) galaxy (see Fig. |3J|. 

3.3 Background galaxies 

The way we simulate the background source population is 
the same for the 2 cases when the lenses belong to a cluster 
versus when they are field galaxies: 

• they are allocated random positions 

• number counts are generated in consonance with galaxy 
counts typical for a 2 hour integration time in the R-band. 
The magnitudes are assigned by drawing the number count 
observed with the CFHT 

• the shapes are assigned by drawing the ellipticity from 
a gaussian distribution similar to the observed CFHT ellip- 
ticity distribution (see Fig. 




dot-dashed: king,846km/s,r core =0. 1 7", 



a=0.3 



r farcsec") 




Figure 3. The reduced shears ( ) for each profile (Panels 
1 and 3) and differences between them expressed in percentages 
(Panel 2 and 4), for L = L*/10 and L = 3L* respectively, where 
once again the PIEMD is the reference profile. 



• Redshift distribution: We use the R-band to define the 
number counts of galaxies and use the HDF prescription in 
terms of the mean redshift per magnitude bin, and the same 
redshift distribution as BBS. 
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magnitude ellipticity 



Figure 4. Galaxy counts from CFH12k and HST data (Left panel), and the ellipticity distributions for the CFH12k and HST data 
(Center panel). The value in parentheses is the mean value of the ellipticity in each case. Right panel: the redshift distribution of the 
background population in our simulations. 



3.4 Lens galaxies 

The way we simulate the foreground lenses galaxies is dif- 
ferent in the case of galaxies inside the cluster and in the 
field. 



3.4-1 Cluster Galaxies 

We put the individual lenses constituting a cluster at a red- 
shift of 0.2, and model it as a superposition of large-scale 
smooth cluster component and a few clumps. In order to 
obtain a better match to the real data, the positions and 
the magnitudes of the foreground cluster galaxies are drawn 
from the positions and the magnitudes of ellipticals in the 
cluster A1689 at z = 0.18. 



3-4-2 Field Galaxies 

For simplicity, the lenses are randomly distributed in posi- 
tion and uniformly distributed in a redshift range from 0.2 
to 0.5. This distribution is a crude approximation of reality. 



4 MAXIMUM LIKELIHOOD ANALYSIS 
4.1 Methodology 

Using the foreground cluster and field galaxies as the lens 
for the sheet of generated background galaxies, we use 
LENSTOOL to solve the lensing equation and produce a 
catalogue of lensed background galaxies. This catalogue con- 
tains the following information for each lensed object: the 
position, the shape parameters, and the redshift. Then this 
catalogue is processed through a numerical code that re- 
trieves the input parameters of the lenses using a max i mum - 
likelihood method as proposed bv lSchneider fc Rbd jl997l) 
and as implemented by Natarajan & Kneib (1997). For each 
image (i), given a mass model for the foreground lensing 
galaxies, we can compute the amplification matrix a» as 
a sum of the contribution from all the foreground galaxies 
j; Zj < Zi that lie within a circle of inner radius Rmin, and 
outer radius R max measured from the centre of the image 
(i): 



E 



(32) 



Zj < Zi 

d(i,j) < R„ 



The total shear experienced by a background galaxy 7i can 
be obtained by summing the contributions from all the fore- 
ground galaxies j;Zj < Zi that lie within an annulus with 
inner and outer radii respectively at Rmin and R m ax- 

Given the observed ellipticity e l obs (defined as e — (a — - 
b)/(a + b)) and the associated amplification matrix at, we 
are able to retrieve the intrinsic ellipticity ef of the source 
before lensing: 



(33) 



In the weak lensing regime, this relation can be simplified 
as: 



+ 7» 



(34) 



In order to assign a likelihood to the parameters used to de- 
scribe the lensing galaxies, we use P" , the ellipticity proba- 
bility distribution in the absence of lensing. Repeating this 
procedure for each image in the catalogue, we construct the 
likelihood function: 



(35) 



which is a function of the parameters used to define the mass 
models of the lenses. For each pair of chosen parameters, we 
can compute a likelihood function. The larger the likelihood 
function, the more accurate the retrieved parameters used to 
describe the lenses. The inversion from the observed elliptic- 
ity to the intrinsic ellipticity is fully analytic and takes into 
account all the non-linearities arising in the strong lensing 
regime, which may occur in the configuration with a cluster 
component. 

The likelihood function for the parametric mass model 
for the lenses does have interesting convergence properties. 
The likelihood surface is topologically complex since the de- 
generacies in the problem tend to produce several maxima. 
However, the convergence in both the cluster lens case and 
field lens case is driven essentially by the width of the intrin- 
sic ellipticity distribution of the sources. The degeneracies in 
this scheme are the generic ones that plague all lensing anal- 
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yses, the mass-sheet degeneracy (essentially the addition of 
a constant sheet of mass to the lens plane does not pro- 
duce any discernable effect on the lensing of the background 
sources), the shot noise due to the finite number of sam- 
pled background sources and the details of the truncation 
of the mass profile of the lenses. The mass sheet degeneracy 
cannot be circumvented since we are necessarily measuring 
both the magnification k and the shear 7 from the same data 
points. Note that these are implicitly needed in computing 
the amplification matrix. Shot noise is more of a limitation in 
ground-based surveys when the number count of background 
galaxies is sparser compared to the space based data despite 
the fact that lensing helps by magnifying fainter sources that 
normally would not make it into a magnitude limited survey. 
The details of the mass profile, and in fact, the prescription 
used to truncate the mass at large radii influences the like- 
lihood results. Since in galaxy-galaxy lensing (both in the 
field and interior to a cluster) we are most sensitive to the 
mass enclosed within an aperture and are less sensitive to 
the density profile in the inner regions the truncation of 
the mass distribution at large radii drives the convergence 
of the likelihood function. This can be clearly seen in the 
plots where the influence of i? max shows up directly in the 
likelihood contours (see Fig. |SJ . Note that the parameters 
chosen to characterise the mass model: the central velocity 
dispersion and the truncation radius also contribute to the 
noise in the likelihood, thereby pointing to more efficient re- 
parameterisations that we also explore in the final section 
of this paper. 



5 RESULTS 

We present the results obtained for the simulated data set, 
for the PIEMD, NFW and truncated power law (PL) mod- 
els for two different configurations fFig. l9l [TU1 and [Til at the 
end of the paper). The points mark the value of the INPUT 
parameters used in order to generate the simulated cata- 
logue, and the cross stands for the value of the OUTPUT 
parameters as estimated from the maximum likelihood anal- 
ysis. We will refer to the different plots by assigning them 
a number, the first one is the upper left plot, and the last 
one (plot number 9) is the lower right plot. The first plot 
shows the reference field situation: 25000 elliptical sources 
in a field of 26x26 arcmin 2 , which translates into a number 
density of about 35 galaxies per arcmin 2 . Then the follow- 
ing plots (2, 3, 4) show the results obtained with 25000 cir- 
cular sources, then with 40000 elliptical sources mimicking 
typical space observations with a density of 60 galaxies per 
arcmin 2 , then with 12500 elliptical sources, corresponding to 
the ground based configuration with a galaxy density of 17 
per arcmin 2 , in a field configuration. Note that in the case 
of circular sources, we do not have to deal with the intrin- 
sic ellipticity noise: the detection is therefore improved and 
the contours are tighter. The plots number 5 and 6 demon- 
strate the effect of the unknown redshift distribution for the 
background sources, in a field configuration. The last row 
represents the cluster configuration: the 'standard' configu- 
ration (plot 7), then configurations where an uncertainty on 
the cluster has been introduced (plot 8 and 9). for each plot, 
the contours represent the 3a,4o",5a confidence levels, and 
along the dotted lines, the mass within a projected radius 
Raper = 100 kpc is constant, equal to the value indicated 
on the plot. 



4.2 Cluster weak lensing mass estimates 

We construct a composite mass model for the cluster by 
superposing a large-scale smooth mass component and indi- 
vidual galaxies. As a first guess for the smooth mass model 
we use the averaged shear field obtained by simply binning 
up the shear in radial bins from the centre outward. This is a 
prior in the analysis that gets modified with every iteration 
once the clumps are added to the model. We simulate the 
deformations induced by a clump with known parameters 
(which can be easily derived from a weak lensing analysis for 
example) . Then we add in the individual cluster galaxies and 
derive the shear of this composite system, which turns out to 
be larger than the shear for the clump alone as expected: this 
implies that we need to simultaneously modify the smooth 
component and the clumps during the optimisation process. 
In massive lensing clusters, about 10% of the total mass i s 
associated with the individual galaxies (Nataraian l2002al) . 
This large scale clump is described by a PIEMD profile with 
the following parameters: cro = 1070 kms , r cu t = 930fcpc 
and r core = 60 kpc. This gives a total integrated mass of 7.3 
10 14 Mq. The mass we partition to galaxies is of the order 
of 7.3 10 13 Mq. How do we modify the large-scale clump's 
parameters in order to match the shear field? We find that 
the velocity dispersion of the large-scale clump needs to re- 
duced by about 5%, keeping the others parameters fixed to 
accomodate the dumpiness. 



5.1 The number of background lenses 

For each profile, in the field configuration, we explore the 
influence of the background density on the detection. The 
'standard' configuration has 25000 background sources (a 
density of 35 galaxies per arcmin 2 ), we explore what happens 
when we increase this number to 40000 (60 per arcmin 2 ), or 
reduce it to 12500 (17 per arcmin 2 ). 17 galaxies per square 
arcminute corresponds to simulating ground based survey 
data, whereas 60 galaxies per square arcminute corresponds 
to the space based survey data. The main difference between 
the ground and space configurations is that from space, the 
statistics are significantly improved, and so the detection 
contours are significantly narrower. 

5.2 The effect of assigning redshifts from an 
assumed distribution 

To quantify the uncertainty arising from not knowing the 
redshifts for background sources, we performed the analy- 
sis after assigning redshifts drawn from a distribution. The 
25000 sources are put at a mean redshift z s and images are 
simulated. When constructing the simulated catalogue, the 
background objects are assigned a mean redshift of z s + S z . 
This catalogue is then input into the maximum likelihood 
code. Since the strength of the shear is proportional to the 
distance between the sources and the lenses, under estimat- 
ing the source redshifts leads to an overestimate of the lens 
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Sz s (%) 


SE 


Sao 


CO 











[174-196] 


20 


0.05 


-7 


[170-190] 


-20 


-0.05 


7 


[180-208] 



Table 1. Influence of an uncertainty in the mean redshift of the 
sources: an error of 8zg correspond to a variation of E(z^,zg) 
equal to 8E, which correspond to a variation on of order 5oq. 
We see from this table that this estimated variation on <ro is 
coherent with the variation in the detection range as derived from 
the maximum likelihood analysis. 



masses and tends to shift the confidence contours toward 
higher values for the velocity dispersion. For the same rea- 
son, systematically overestimating the source redshifts leads 
to an underestimate of the galaxy masses. In any case, we 
find that a redshift uncertainty of ±0.2 does not dramati- 
cally modify the conclusions. This is the typically the preci- 
sion we can get with photometric redshift estimation which 
is encouraging for future surveys. This results is coherent 
with a recent study by Kleinheinrich et al., 2004 based on 
the galaxy-galaxy lensing results from the COMBO-17 sur- 
vey. They found that it is of great importance to know the 
redshifts of individual lens galaxies in order to constrain the 
properties of their dark matter halos, but that the knowl- 
edge of individual source redshifts improves the measure- 
ments only very little over use of statistical source redshift 
distribution. 

To be more quantitative, let us consider the lensing 
equation and express it for a constant deflection angle. For 
a PIEMD profile with a given ao, we have: 

a o —pr S — constant; and J^ S = E{zi,Z s ) (36) 
D.s Ds 

as introduced in IGolsd (2002). Therefore the equation can 
be rewritten as: 

(JqE(zl, zs) = constant (37) 

The lenses are kept at a redshift of 0.2, and the mean red- 
shift of the sources is changed by 8z = ± 20 %. We then 
evaluate the corresponding Sao error introduced in the re- 
trieval of the central velocity dispersion in the likelihood 
analysis. Table Ogives the results: when we put an error of 
Szs (%), this give a variation of E(zl, zs) equal to SE, and 
the corresponding variation on ao is of order Sao- This range 
of values is given by the projection of the 3a contours along 
the ao axis. We can see that the variation in the detection 
range is in agreement with the calculations made. 

5.3 Influence of the uncertainty in the cluster 
modeling 

When working with the real data, we will have to put in 
by hand the description of the cluster. The reliability of 
the results will depend on the accuracy with which we de- 
scribe cluster. In order to study the influence of the uncer- 
tainty of the cluster profile, the cluster component is de- 
scribed by a PIEMD profile with a velocity dispersion of 
a duster = 1000 kms -1 . When constructing the simulated 



catalogue, the cluster component is assigned a velocity dis- 
persion of a c i us ter + Str. The likelihood is then computed for 
this case. 



6 DISCUSSION AND CONCLUSIONS 
6.1 Constraints obtained on mass profiles 

For the PIEMD and NFW profiles, we found that we are 
able to retrieve the characteristic halo parameters with ex- 
tremely good accuracy, for every configuration. In fact, in- 
terestingly enough, the dotted lines in Figs. 9, 10, and 11 
show us that the aperture mass is retrieved very accurately. 
This immediately suggests the re-parameterisation of mod- 
els considered here: rather than fitting in the (ao,r) plane, 
we can fit the deformations directly in the (M aper , R ap er) 
plane. This formulation is explored in the next section. 

For the PL profile, we find that we can put some con- 
straints on a, the slope of the density profile, but not on 
the velocity dispersion, since the likelihood function does 
not always converge along that direction. So we can use this 
profile to estimate the slope of dark matter halos, without 
trusting the constraints we get on the velocity dispersion for 
the profile. 



6.2 Influence of Rmax 

When working with real data, the results we get depend on 
the value chosen for R ma x- When we take a low value for 
Rmax, the shear for the image (i) is calculated with fewer 
lenses, and we find that the contours do not close in r cut . 
When the value of this parameter is increased , the contours 
converge and close. On the other hand, picking a high value 
for Rmax introduces some noise in calculating the shear, 
and can dilute the galaxy-galaxy lensing signal significantly: 
lenses that do not effectively participate in the lensing of an 
image if utilised in the calculation become a source of noise. 
Fig-EJillustrates precisely this situation: for a given value of 
Rmax, we obtain a good estimate of ao, but the robustness 
of the constraint on r cu t is directly related to the value of 
Rmax- Others authors have reported that their results are 
sensitive to the value of this parameter (e.g. Kleinheinrich, 
PhD thesis). 

The choice of Rmax is therefore important. To get the 
order of magnitude of this parameter, we compare the char- 
acteristic noise in the problem, i.e. ^7=, to the signal we are 
sensitive to i.e. the reduced shear - the factor of 0.25 is the 
width of the intrinsic ellipticity distribution and N is the 
number of background objects at a distance r from a lens- 
ing galaxy. This noise has been estimated by analysing data 
from the ground-based CFH12k observations of the cluster 
A1763 at a redshift of z = 0.22. From space, we expect 
the number of background objects to be about 6/7 times 
higher. Fig. HJ shows that a value of about 100" can be used 
for R 

ma x • 

On the other hand, the choice of the parameter R m in 
does not influence the results, so we fix this parameter to be 
of the order of a few kpc. 
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Figure 5. The PIEMD profile, in a field configuration: In the 
panels from left to right, the value of R m ax increases from 30", 
60" , 90" to 150" . When R max > 100" , the contours coverge along 
the r cu t axis. Note that the convention throughout this work is 
that the dot marks the value of the input parameters and the 
cross marks the retrieved output values. 
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Figure 6. Comparing the signal to noise for an estimate of the 
optimal choice for Rmax - The solid lines correspond to the charac- 
teristic cllipticity (due to the width of the ellipticity distribution), 
as a function of radius, from ground and space. The other 
line types (dot, dashed and dot-dashed) correspond to the signal - 
the reduced shear as a function of radius for typical luminosities. 



likelihood, for instance M ap er and R ap er, which we explore 
below. 

We have: 

C — C(ao,r cut ) & M ap er = Maper {Raper, &0,r cut) (38) 

so we can write: 

O"0 = 0"o(M aper , Raper, Tcut) (39) 

the likelihood function then becomes: 

f~-{^Iaper , Raper , i" cut) (40) 

and by summing over r cu t, we get: 

C! = 2_j C-{Maper, Raper, T cut) (41) 
r C ut 

C — C (A'laper , Raper) (42) 

The sum is performed for a set of r cu t values around 
the input value used to simulate the catalogue and the range 
denned by the projection of the 3a contour level along the 
r C ut axis. The results do not depend strongly on the range 
used to do the sum. Fig. |7| and [8] shows the £.' contours 
we get for the NFW and the PIEMD profiles. Since the PL 
profile does not have a cut-off radius we cannot compare it 
easily with the re-parameterized PIEMD and NFW profiles. 
The plots above show that we can put strong constraints 
on the aperture mass, the crossed line represents the line 
M a per{Raper) as computed with the input model used to 
generate the simulated catalogue. 

The motivation of such a re-parameterisation is that we 
deal with more direct physical quantities than halo parame- 
ters, that is an aperture mass calculated within an aperture 
radius. The primary motivation for galaxy-galaxy lensing 
studies was to measure halo masses, so this offers a more 
convenient parameterisation for that purpose. This is also a 
different way of measuring masses compared to the aperture 
densitometry method. Moreover, it is more suited to the case 
of clustered galaxies since we are not able to integrate the 
shear profile for any individual galaxy. 

To conclude in this paper, we have discussed galaxy- 
galaxy lensing in the context of measuring masses of field 
and cluster galaxies. We compare the robustness of recov- 
ering input parameters for the mass distribution of lenses 
from a ground based survey and space based observations. 
We explore a wide range of input mass models for galaxy 
halos. We simulate the galaxy-galaxy lensing effect and gen- 
erate synthetic catalogues. A maximum likelihood method is 
applied to the catalogues to successfully recover the lens pa- 
rameters in various configurations. Going beyond the stan- 
dard parameterisation of a dark matter halo, we propose a 
re-parameterisation of the problem in terms of more direct 
physical quantities: the aperture mass calculated within an 
aperture radius. The main result of this re-parameterisation 
is that we are able to put even stronger constraints on the 
aperture mass for an L* galaxy. 



7 RE-PARAMETERIZATION OF THE 
PROBLEM 

Thus far, we have performed the likelihood analysis to op- 
timize the values of two parameters a and r cu t (or r s ). A 
different set of parameters can be chosen for maximizing the 
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Figure 7. The likelihood C'(M aper , R ap er) for the NFW (up- 
per panel) and the PIEMD profile (lower panel). The crossed line 
represents the M aper (R aper ) contour obtained with the input pa- 
rameters used to describe the foreground lenses in the cluster 
configuration. 
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Figure 8. The likelihood C (M aper , R aper ) for the NFW (up- 
per panel) and the PIEMD profile (lower panel). The crossed 
line represents the M aper (R aper ) contour obtained with the in- 
put parameters used to describe the foreground lenses in the field 
configuration. 
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Figure 9. Results for the PIEMD profile. The top left plot (number 1) is the reference field situation with 25,000 elliptical background 
sources in a field of 26x26 arcmin 2 , corresponding to a density of background sources equal to about 35 galaxies per arcmin 2 . Then the 
following plots (2, 3, 4) show the results obtained with 25,000 circular sources, then with 40,000 elliptical sources mimicking typical space 
observations with roughly 60 galaxies per arcmin 2 , then with 12,500 elliptical sources, corresponding to the ground based configuration 
with a galaxy number density of 17 galaxies per arcmin 2 . The plots number 5 and 6 show the effect of introducing an uncertainty in the 
mean redshift of the source population: this uncertainty is equal to -20% (left) and +20% (right). The last row show the results from 
the cluster configuration: the reference cluster configuration (plot 7), then configurations where an uncertainty on the cluster modeling 
is introduced: knowing the central velocity dispersion of the cluster to within ±10%. The contours in this figure represent the 3cr,4a,5cr 
confidence levels, and along the dotted lines in every panel, the mass within a projected radius R ape r of 100 kpc is kept constant at the 
value indicated on the plot. Note that the dot indicates input values and the cross the retrieved output. 
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Figure 10. Results from the NFW profile: the top left plot (number 1) is the reference field situation with 25,000 elliptical background 
sources in a field of 26x26 arcmin 2 , corresponding to a density of background sources equal to about 35 galaxies per arcmin 2 . The 
following plots (2, 3, 4) show the results obtained with 25,000 circular sources, then with 40,000 elliptical sources, corresponding to data 
obtained from space with 60 galaxies per arcmin 2 , then with 12,500 elliptical sources, corresponding to the ground-based data of about 
17 galaxies per arcmin 2 . The plots number 5 and 6 show the effect of introducing an uncertainty in the mean redshift of the source 
population: this uncertainty is equal to -20% (left) and +20% (right). The last row show the results from the cluster configuration: the 
reference cluster configuration (plot 7), the effect of introducing an uncertainty of ±10% in the central velocity dispersion of the cluster 
model used. Contours in these figures represent the 3<r, 4cr, and 5<r confidence levels, and along the dotted lines in each panel, the mass 
within a projected radius R aper of 100 kpc is kept constant, with M apor at the quoted value indicated on the plot. Note that the dot 
indicates input values and the cross the retrieved output. 
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Figure 11. Results from the PL profile. The top left plot (number 1) is the reference field situation with 25,000 elliptical background 
sources in a field of 26x26 arcmin 2 , corresponding to a density of background sources equal to about 35 galaxies per arcmin 2 . Then the 
following plots (2, 3, 4) show the results obtained with 25,000 circular sources, then with 40,000 elliptical sources corresponding to a 
typical space based observation yielding 60 galaxies per arcmin 2 , then with 12,500 elliptical sources corresponding to ground based data 
with a galaxy number density of about 17 galaxies per arcmin 2 . The plot number 5 and 6 show the effect of introducing an uncertainty in 
the mean redshift of the source population: this uncertainty is equal to -20% (left) and +20% (right). The last row show the results from 
the cluster configuration: the reference cluster configuration (plot 7), and then configurations where uncertainty has been introduced in 
the cluster modeling. The contours in all these panels represent the 3cr, 4a and 5<r confidence levels. Note that the dot indicates input 
values and the cross the retrieved output. 



